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Abstract 



Local expansion exponents for nonequilibrium dynamical systems, described by 
partial differential equations, are introduced. These exponents show whether the 
system phase volume expands, contracts, or is conserved in time. The ways of calcu- 
lating the exponents are discussed. The principle of minimal expansion provides the 
basis for treating the problem of pattern selection. The exponents are also defined 
for stochastic dynamical systems. The analysis of the expansion-exponent behaviour 
for quasi-isolated systems results in the formulation of two other principles: The 
principle of asymptotic expansion tells that the phase volumes of quasi-isolated sys- 
tems expand at asymptotically large times. The principle of time irreversibility 
follows from the asymptotic phase expansion, since the direction of time arrow can 
be defined by the asymptotic expansion of phase volume. 
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1 Introduction 



Macroscopic properties of nonequilibrium dynamical systems are usually characterized 
by the Gibbs entropy defined for the appropriate probability measure [1,2]. The usage 
of this entropy has two shortcomings: First, it is defined only for steady states, and 
an effective definition of entropy far from equilibrium is not known [2]. Second, there 
appears the problem of singularity resulting from the following [1,2]. A steady state of 
a dynamical system corresponds to the motion on an attractor. The related probability 
measure is, typically, singular, because of which the entropy becomes equal to — oo. To 
avoid this problem, one considers not the entropy itself but the entropy production or 
entropy production rate, such as the Kolmogorov-Sinai entropy rate [3]. The latter are 
finite, but again, they are defined only for stationary states. 

Employing the Gibbs entropy in nonequilibrium statistical mechanics confronts as 
well another paradox. This entropy, for an isolated system, as is well known, does not 
change in time. While the second law of thermodynamics requires that the entropy of an 
isolated nonequilibrium system would be an increasing function of time. To overcome this 
contradiction, one introduces modified variants of entropy, among which the most popular 
is a coarse-grained enropy [4]. Other variants are the von Neumann entropy [5] and the 
average entropy [6], which change with time for an isolated system not in equilibrium. 
Lebowitz [7] argues that the relevant entropy for understanding the time evolution of 
macroscopic systems is not the Gibbs entropy Sq = — Trp In p but the Boltzmann entropy 
Sb = In |T|, where T is the system phase volume. Contrary to the Gibbs entropy, which 
is constant for an isolated macroscopic system, the Boltzmann entropy should typically 
increase with time. A weak point in resorting to the Boltzmann entropy Sb is that it is 
not always clear how to determine the phase volume for a strongly nonequilibrium system. 

Instead of trying to construct a special form of entropy that would be valid for nonequi- 
librium dynamical systems, it is possible to employ another characteristic. The aim of the 
present paper is to demonstrate that the local expansion exponent can be such a convenient 
quantity characterizing the temporal evolution of dynamical systems. In Section 2, the 
local expansion exponent is introduced for an infinite-dimensional dynamical system and 
its properties are described. Section 3 shows the usefulness of the expansion exponents 
for treating the problem of pattern selection. The generalization of the expansion expo- 
nents to stochastic dynamical systems is given in Section 4. Several examples of Section 
5 illustrate how the exponents can be calculated. The notion of quasi-isolated systems 
is formulated in Section 6 and some particular cases are considered. Section 7 contains 
conclusions and discussion. 

2 Local expansion exponents 

Let x G D C M. d be a (/-dimensional set of continuous variables given on a domain D and 
let t G IR+ = [0, oo) denote time. Suppose we are interested in the evolution of a family of 
complex-valued functions y~i(x,t), with i G N + = {1,2, . . .}, describing a nonequilibrium 
physical process. It is convenient to employ the matrix notation [8], where the pair {i,x} 
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is treated on equal footing as a point in the label space N+ x D. Then the dynamic state 
y(t) = [yi(x,t)\ is presented as a vector-column with respect to the multi-index {i,x}. 
The components yi(x,t) may correspond to some physical quantities. In particular, these 
could be statistical or quantum averages of operators for local observables. This could 
also be a wave function or a density matrix. 

Dynamic state y(t) G T pertains to a phase space T . Let the velocity field v(y,t) = 
[vi(x,y,t)\ be given as a column with respect to the label pair {i,x} and defined in a 
space tangent to T . The set of evolution equations can be compactly written as 

j t y(t)=v(y,t). (1) 

For the dynamical components yi(x,t) G Ti{x) pertaining to the phase subspaces ^(x), 
the structure of the total phase space T is that of the tensor product 

The continuous product of spaces Ti{x), with respect to the continuous variable x, is a 
direct generalization [9] of the discrete product over the discrete index % = 1,2, . . .. The 
related continuous product of functions f(x) is naturally defined [9] as 



II f( x ) = ex P / ln f( x ) dx ' 



with the integration over D, where any branch of the logarithm may be taken. The 
elementary phase volume at time t is interpreted [10,11] as 

sm^nnsy^t). (2) 

i x 

Introduce the local expansion exponent 



X{t) = In 



ST(t) 



*T(0) 



(3) 



which shows whether the phase volume, at time t, expands, contracts, or is conserved. 
Rewriting the definition (3) in the form 

1^)1 = 1^(0)1^, (4) 

one sees that positive X(t) implies expansion, negative X(t) means contraction, while zero 
X(t) tells that the phase volume at time t does not change. Note that X(t) could equally 
be called the local contraction exponent. However, to term it the expansion exponent 
looks more logical as far as X(t) > implies expansion. 

One may notice that the quantity (3) resembles the entropy production S(t) — S(0), 
with the entropy S(t) = \n\5T(t)\ defining a la Boltzmann. But, clearly, this is just a 
resemblance since the expansion exponent (3) has sense for any dynamical system and 
does not require to be ascribed to a statistical system. For the latter, one usually defines 
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the entropy production being based on the Gibbs entropy, which allows for the description 
of only steady states [1,2]. But for strongly nonequilibrium systems the Boltzmann form 
is more appropriate [7]. 

The expansion exponent (3) is directly related to the multiplier matrix M(t) = 
[Mij(x,x',t)] with the elements 

and with the initial condition 

M i:j (x, x', 0) = Sij S(x - x') . (6) 
Invoking equations (2) and (5), we have 

X (t) =E f^\Mii(x,x,t)\ dx. (7) 
i J 

This expression can be rewritten in a general form not depending on a particular repre- 
sentation. For this purpose, let us recall that for a matrix A the equality 

Tr In A = In deti (8) 

is valid. The proof of this equality is based on the fact that the trace of a matrix does 
not depend on a matrix representation. Then, employing the diagonal representation for 
A = [5 mn A n ], one has 

Tr In A = ln A n = In Y[ A n = In deti , 

n n 

which proves equality (8). Defining the matrix 

ln|M(t)| = [ln| M^x.x'M , 
equation (7) may be presented as 

X(t) = Trln | M(t) | . (9) 

And using equality (8), we obtain 

X(t) =ln|detM(t)| . (10) 

Thus, the expansion exponent (3) is expressed through the multiplier matrix (5). 

Another expression for the expansion exponent can be found by invoking the evolution 
equation (1). Define the Jacobian matrix J(t) = [Jij(x,x',t)} with the elements 

t ( ' +\ - 8vi(x,y,t) , 



5 



The variation of equation (1) over the initial conditions results in the equation 

j t M(t) = J(t)M(t) (12) 

for the multiplier matrix. The initial condition for equation (12), according to equality 
(6), is 

M(0) = i = [5ij S(x - x')} . (13) 
With the usage of equation (12), we have 

4- TrlnM(t) = TrM _1 (t) ^ M(t) = TrikT 1 ^) J{t)M(t) = TrJ(t) , 
at at 

from where it follows that 

4- Tr In M(t) = Tr J(t) . 
at 

Integrating this equation gives 

TrlnM(t) = f Tr J it') dt' . 
Jo 

And equality (8) yields 

TrlnM(t) = lndetM(t) . 
Comparing the latter two equations results in the equality 

detM(t) = exp j jf* Tr J(t') dt'} . (14) 

Substituting expressions (14) into equation (10) provides us with the local expansion 
exponent 

X(t) = Re f Tr J{t') dt' . (15) 



Equations (10) and (15) are the main forms for calculating the expansion exponent. 

For practical calculational purposes, it is possible to use different representations con- 
nected with a chosen basis {(p n (t)} of the vector-columns (p n (t) = [ip n i(x, t)]. In what 
follows, such a basis is assumed to be orthonormal and complete, 

¥>m(*)Vn(f) = Srnn , ^ ip n (t)lfi+(t) = 1 . 

n 

From here, we get a useful equality 

^fVw + rfi^-O (16) 

for any indices m and n. The n-representations for the multiplier and Jacobian matrices 
are defined by the elements 

M mn {t) EE <p+(t)M(t)<f n (t) , J m n(t) = ^(t)J(t)<f n (t) . (17) 
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These elements are connected with each other through equation (12) yielding 

' , (18) 



It M -" W = ? 



J mk (t)M kn (t) + M mk (t)<pt(t) djf ^ - M kn (t) 



where equality (16) has been invoked. The initial condition for equation (18) is M mn (0) = 
5 mn . The n-representation can be exploited for deriving several important properties. 

Proposition 1. If the multiplier matrix M(t) possesses the eigenvalues (J> n (t), given 
by the eigenproblem 

M(t)tp n (t) = fi n (t)<p n (t) , (19) 
with the eigenvectors composing a complete orthonormal basis, then the eigenvalues 

/i n (t) = exp j^J n „(t') dt'j (20) 

are expressed through the diagonal elements J nn (t) of the Jacobian matrix. 

Proof is given in Appendix A. Note that the Jacobian matrix is, in general, nondiag- 
onal. 

Proposition 2. Suppose a complete orthonormal basis {<p n (t)} is such that 

^ = (m + n) . (21) 

Then <p n (t) are the eigenvectors of the multiplier matrix M(t) if and only if they are also 
eigenvectors of the Jacobian matrix J(t), with the corresponding eigenvalues fJ, n (t) and 
J n {t) connected by the relation 

Ai n (*)=exp{jT* J n (0*'} . (22) 

Proof is in Appendix B. The difference between equations (20) and (22) is in the 
form of the Jacobian matrix. In equation (20), J n n{t) is a diagonal element of a generally 
nondiagonal matrix. While in equation (22), J n (t) is an eigenvalue of the Jacobian matrix, 
which hence is diagonal. In other words, theorem 2 tells that, if condition (21) is satisfied, 
then the multiplier matrix is diagonal if and only if the Jacobian matrix is diagonal. 

Proposition 3. Let the Jacobian matrix have the limit 

lim J(t) = J (23) 

and let J n be the eigenvalues of J. Then the limit 

Ito A. (24) 
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is the sum of the Lyapunov exponents 



A n = ReJ„ = lim - In \n n (t)\ . 



(25) 



Proof is in Appendix C. The sum of the Lyapunov exponents gives the entropy pro- 
duction rate for stationary states [1,2]. Hence, the left-hand side of equation (24) has the 
meaning of the entropy production rate. Recall that summing up only positive Lyapunov 
exponents leads to the Kolmogorov-Sinai entropy rate [3]. Contrary to the entropy pro- 
duction, defined on the basis of the Gibbs entropy [1,2], the local expansion exponent (3) 
does not meet the problem of singular measure and it has sense not only for stationary 
states but for any strongly nonequilibrium states. 

3 Problem of pattern selection 

The local expansion exponent turns out to be a crucial notion for solving the problem of 
pattern selection. This problem arises when the evolution equations for a physical process 
possess a multiplicity of solutions corresponding to different spatio-temporal structures 
[12]. When all these solutions are stable, it is not clear how one could distinguish between 
them and, respectively, what spatio-temporal structures should be treated as more prefer- 
able. At the same time, nature does prefer some structures over the others, since some 
of them arise more often for the given experimental protocol. This problem of pattern 
selection can be solved as follows. 

Let the multiplicity of solutions be parametrized by a label (3 pertaining to a label 
manifold B = {{3}. The submanifodl of B related to stable solutions is called [12] stability 
balloon. Each dynamic state y(/3,t), labelled by (3, corresponds to a particular spatio- 
temporal pattern. Because of the parametrization of dynamic states by the label /3, the 
dependence on this label also comes through the multiplier matrix M((3,t), Jacobian 
matrix J(/3,t), and the local expansion exponent X(/3,t). Since nature distinguishes 
among possible patterns differentiating them on more or less preferable, there should 
exist a probability distribution classifying these patterns as more or less preferable. Such 
a probability distribution of patterns p((3, t) can be obtained by minimizing an information 
functional 



under the normalization condition on the probability distribution, where the integration 
is over the label manifold B. 

The information density I p (/3,t) can be found from the following requirements: The 
existence of an invariant probability measure 




p(f3,t)\5T((3,t)\=p([3,0)\5T(f3,0)\ 



(26) 



and the occurence of an invariant information measure 



i p {p,t)\sr{i3,t)\ = i p (p,o)\6r(p,o)\ . 



(27) 
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These equalities tell us that the probability of a pattern in an elementary phase volume 
and the amount of information in such a phase volume are temporal invariants. We also 
require that at the initial time t — the information density be given by the standard 
Shannon form 

I p ((3,0)=p((3,0)\np((3,0). (28) 
Then, from the invariance of the probability measure (26), we have 



p(P,0)=p(l3,t) 



sr(p,t) 



sr(p,o) 



Substituting this into the Shanon information density (28) and using the invariance of 
the information masure (27), we obtain the information density of patterns 

W t) = p(P, t) [lnp(/3, t) + X(P, t)\ . (29) 

Thus, the information functional to be minimized is the pattern information 

I p (t) = J p((3, t) lnp((3, t)df3 + J p((3, t)X((3, t) df3 . (30) 

The minimization of the pattern information (30), under the normalization condition on 
p(/3,t), implies the minimization of the conditional information 



I p (t)=I p (t)+l(t) \jp(p,t)d(3- 



(31) 



in which 

l(t) = \nZ(t) - 1 

is a Lagrange multiplier. Minimizing information (31) results in the pattern distribution 

pi(3 ' t) = W) exp{-x(/M)} ' (32) 

with the normalization factor 

Z(t) = J exp{-X((3,t)}d(3 . 

The expansion exponent, according to equations (10) and (15), can be determined by one 
of the forms 

X(/3, t) = In |detM(/3, t) \ , X(f3, t) = Re /' TrJ(/3, t') dt' . (33) 

Jo 

With the first of these forms, the pattern distribution (32) becomes 

p(P,t)= , (34) 

V ; \detM(P,t)\ V ; 
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where the normalization factor is 



Z(t)=/ 



dp 

\detM((3,t)\ 



The expression for the pattern distribution (32) shows that those patterns are more prefer- 
able whose expansion exponents are smaller. The most probable pattern corresponds to 
the minimal expansion exponent, 

max p(/3,t) mmX(P,t) , (35) 

which constitutes the principle of minimal expansion. 

The probabilistic approach to pattern selection was advanced in Refs. [10,11], where 
the pattern distribution was rather postulated by invoking analogies with statistical me- 
chanics. Here we have shown how this pattern distribution can be derived from the 
minimization of the pattern information. The probabilistic approach to pattern selection 
was applied to the problem of turbulent photon filamentation [11,13], with theoretical 
description being in a very good agreement with all experimental observations. This ap- 
proach makes also the basis for the method of self-similar prediction [14-16] applied for 
analysing and forecasting possible scenarios in the behaviour of nonequilibrium complex 
systems, such as markets. In the latter case, the scenario distribution of form (34) was 
employed. 



4 Stochastic dynamical systems 

The notion of expansion exponents can be generalized to stochastic dynamical systems. 
Suppose the evolution equations include a stochastic field = [£i(x, t)]. Then the 
dependence on this field enters the dynamic state y(£,t) = [yi(x,£,t)] as well as the 
velocity field v(y,£,t) = [vi(x,y,£,t)]. All these quantities are again treated as columns 
with respect to the multi- index {i, x} spanning the label space N + x D. The set of evolution 
equations acquires the form 

j t y(Z,t) = v(y,Z,t), (36) 

with the initial condition 

l/(e,0)=y(0). (37) 

The observable quantities are obtained after averaging over the stochastic fields. Denoting 
the stochastic averaging through the double angle brackets <C . . . 3>, we have 

y(*) = «y(£,*)» . (38) 

In the elementary phase volume (2), we now need to set 

8yi(x, t) = « 8yi(x, £, f) > . (39) 

Then the same definition (3) for the local expansion exponent holds true. 
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Introducing the stochastic multiplier matrix M(£,t) = [Mij(x, x', £, t)], with the ele- 
ments 

U «^'M'W$ (40) 

and the initial condition 

M ij (x,x',Z,0)=5 ij 8(x-x') , (41) 
we may define the average multiplier matrix 

M(*) = «M(f,t) > . (42) 

Then expressions (7) and (10) for the expansion exponent retain their sense, with the 
multiplier matrix (42). 

The stochastic Jacobian matrix is defined as J(£,t) = [Jij(x,x',£,t)], whose elements 

are 

_ 5vj(x,y,£,t) 

The variation of the evolution equation (36) yields the equation 



±M(t,t) = J(t,t)M(t,t) (44) 

for the stochastic multiplier matrix (40), the initial condition being 

M(£,0) = 1. (45) 

The averaging of equation (44) gives 

^M(t)=« J(e,t)M(e,t)» . 

Since the right-hand side here, in general, cannot be factorized, there is no presentation 
for the expansion exponent similar to equation (15). Therefore, for stochastic dynamical 
systems, the form (10) of the expansion exponent remains the main, where the multiplier 
matrix is given by equation (42). 

With the help of an orthonormal complete basis {ip n (t)}, we may pass to the n- 
representation 

M mn (£,t) = <ptM(U)<Pn(t) > Jmn{U) = <pt(t)J(U)<Pn(t) > ( 46 ) 

which is analogous to definition (17). The elements (46) are connected by the same 
equation (18). If the eigenproblem 

M(Z,t)<p n (t) = t i n (Z,t)<p n (t) (47) 

holds true, then the average multiplier matrix (42) possesses the same eigenvectors (p n (t), 
with the eigenvalues 

H n {t) = </i n (£,t) > . (48) 
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It is straightforward to reformulate the theorems 1 and 2 for the stochastic multiplier 
matrix M(£,t): 

Proposition 1. If the stochastic multiplier matrix M(£,i) possesses the eigenval- 
ues fi n (!;,t), given by the eigenproblem (47), with the eigenvectors forming a complete 
orthonormal basis, then the eigenvalues 

/i n (£, t) = exp | jf* J nn (£, f) dt'} (49) 

are expressed through the diagonal elements J nn {€,t) of the stochastic Jacobian matrix. 

Proposition 2. Assume a complete orthonormal basis {<p n {t)} is such that con- 
dition (21) is valid. Then the stochastic multiplier matrix M(£,i) is diagonal in the 
n-representation if and only if the stochastic Jacobian matrix is diagonal, 

their eigenvalues being connected by the relation (49). 

But theorem 3 is not valid for stochastic dynamical systems, since J(£,t), in general, 
is not defined for t — > oo. Thus, to find the expansion exponent, we have to resort to 
formula (10). When the stochastic multiplier matrix satisfies the eigenproblem (47), the 
expansion exponent (10) reduces to the expression 

X(t)= ]T In | «/i n (£, t)»|. (50) 

n 

This is the principal form of the expansion exponent for stochastic dynamical systems. 



5 Exponents for stochastic systems 

Now we shall illustrate how the expansion exponents can be calculated for several examples 
of stochastic dynamical systems. We shall consider the stochastic field £(t) as a real- valued 
random Gaussian variable, which is centered at zero, 

< Z(t) > = o . 

For the case of coloured noise, the correlation function 

« mm » = c(t - o 

can be written [17] as 

C(t) = C exp\ 



\t\ 

T 



Setting Co = (71/r + 27I), we may consider two opposite limits, that of a very short 
correlation time r and of a long correlation time. In the short correlation time limit 
r — > 0, taking into account that 

Y T ex P (- 7) ^ ^) 0- - 0) , 
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we have 

limC(t) = 2 7l 5(t) . 

Such a delta-correlated noise is called white because of the uniformity of its spectral 
function. In the opposite limit of a long correlation time r — > oo, we get 

hmC(t) = 2 72 2 . 

This kind of noise can be called infrared, since its spectral function is centered at zero 
energy. In the general case, one should deal with the coloured noise with a finite correlation 
time. To simplify calculations, at the same time keeping tracks of two admissible limits 
of short and long correlation times, one can accept the correlation function of the form 

C(t) = 2 7l 5(0 + 2 72 2 , 

representing a mixture of white and infrared noises. In what follows, we shall consider 
the Gaussian stochastic fields with the correlation function 

« mm » = 2 7l *(f - o + 2 72 2 , (51) 

which includes both limiting cases. Setting here 7 i — > 0, we get the infrared noise, while 
for 7 2 — > 0, we come to the white noise. 

Stochastic differential equations, as is known, can be interpreted either in the sense of 
Ito or in the sense of Stratonovich [17,18]. We prefer to deal with the Stratonovich inter- 
pretation, which is better motivated physically. Another possibility could be to employ 
the stochastic expansion technique [19,20], by presenting stochastic fields as expansions 
over smooth functions with random coefficients. This method allows for the usage of the 
standard differential and integration analysis. The final results of the expansion tech- 
nique coincide with the corresponding expressions derived by means of the Stratonovich 
method. 



5.1 Oscillator with random attenuation 

Let us start with a case of an one-dimensional dynamical system described by the evolution 
equation 

^L = (iu-r + t)y, (52) 

in which y — y(£,,t), £ = £(£), iv and T are positive parameters. This equation corresponds 
to an oscillator with random attenuation. The Jacobian related to equation (52) is 

J(t,t)=w-T + t(t), (53) 
because of which the multiplier is 

, t) = exp {(iw - T) t + J* £(f) dt'} . (54) 
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For the Gaussian random variables, one has the averaging property 

<exp|a^ < e(t')^ , }> = exp|y^<e(0e(O> dt' tft"j , 

where a is a complex number. Then for the local expansion exponent (50), we find 

X(t) = -rt + 7l i + 7 !t 2 . (55) 

If there is no noise, that is 7 i = 7 2 = 0, the motion is contracting. But in the presence of 
noise such that either 7i > T and 7 2 = or 72 7^ for any 7 i, the dynamics is expanding. 
As is seen, the existence of noise can drastically change the system dynamics. 

5.2 Stochastic diffusion equation 

Consider an infinite-dimensional dynamical system presented by the diffusion equation 
with a random diffusion coefficient, 

§H^)§. («) 

where y = y(x,£,t), £ = £(t), and D > 0. For the spatial variable :r given on a finite 
interval, it is always possible, by an appropriate scaling, to reduce this interval to the unit 
one, so that x G [0, 1]. The initial condition is 

y(x,£,0)=y(x,0). (57) 

And let the boundary conditions be 

y(0,S,t) = bo, y(l,£,t) = b 1 . (58) 
From here, the boundary conditions for the multiplier matrix are 

M(0,x\^t) = M(l,x'^,t)=0. (59) 

Hence, if the multiplier matrix possesses the eigenvectors p n = [<p n (x)], their components 
<p n (x) should satisfy the boundary conditions 

<p n (0) = <p n (l) = . (60) 

For the Jacobian matrix, associated with equation (56), we have 

J( x ^'^) = (D + i) — 5{x-x'). (61) 

The eigenproblem for the Jacobian matrix reads 

fi 

J(x,x',g)<p n (x') dx' = J n (£)<p n (x) . 
14 



The eigenfunctions, satisfying the boundary conditions (60), are 

(Pn(x) = V% sin(7rna;) (n = 1, 2, . . .) . 

And the eigenvalues of the Jacobian matrix are 

J n (0 = -(D + S)* 3 n 2 . (62) 

The eigenvectors </?„, not depending on time, satisfy condition (21). Then, by theorem 2, 
the multiplier matrix possesses the same eigenvectors, with the eigenvalues (49), which 
yields 

//„(£, t) = exp S^-n 2 n 2 Dt - vrV J f (*') rft'j . (63) 
In this way, the expansion exponent (50) becomes 



X ( t ) = E [-K 2 n 2 Dt + vr 4 n 4 ( 7l i + 7 f t 2 )] . (64) 



Using the sums 



N AT jV 1U 

sr n i = !l(2N 2 + 3N + l) , £n 4 =^(6Ar 4 + 15AT 3 + 107V 2 -l) , 
we obtain 

X(t) ~ - y JV 3 !* + y iV 5 ( 7l f + 7 2 2 t 2 ) , (65) 

where JV — > oo. In the absence of noise, when 71 = 72 = 0, the system is contracting. 
But if either 7! or 72 are not zero, the phase volume expands. 

5.3 Stochastic Schrodinger equation 

The Schrodinger equation, generated by a Hamiltonian H, with added stochastic back- 
ground, can be written in the form 

= (-M + M , (66) 

which is called normal for dynamical systems. Here ip = i[)(r,£,t), £ = £(£), H = H(r), 
and H — 1. 

For the Jacobian matrix, we have 

J(r,r',0 = H# + 5(r-r') . (67) 
Let the stationary Schrodinger equation be 
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The Jacobian matrix J(£), with the elements (67), possesses the eigenvectors ip = [ip n (x)} 
and the eigenvalues 

J n (S) = -iE n + £. (68) 

Since the eigenvectors ip do not depend on time, condition (21) is valid. Then, by theorem 
2, we find the eigenvalues 

//„(£,*) = exp{-i£ n * + jf*£(f) dt'} (69) 

of the multiplier matrix. And for the expansion exponent (50), we find 

X(t) = ( ll t + 1 2 2 t 2 )N . (70) 

When the noise is absent, the phase volume is conserved. But as soon as either 7! or 72 
are not zero, the dynamical system expands. 



6 Exponents for quasi-isolated systems 

All real physical systems are never completely isolated, but are subject to weak uncontrol- 
lable random perturbations from surrounding. This fact has been repeatedly emphasized 
and discussed in literature [6,9,22-25]. Even more, as has been stressed [26,27], the notion 
of an isolated system as such is logically self-contradictory. This is because to practically 
realize the isolation one has to use technical devices acting on the system, and to ensure 
that the latter is kept isolated, one must apply measuring instruments perturbing the 
system. The preparation and registration processes may essentially influence the system 
evolution [5,28]. The impossibility of isolating macroscopic systems from their environ- 
ments is often considered as the cause of the encrease of entropy required by the second 
law of thermodynamics and the related irreversibility of time arrow [28]. 

To analyze the influence of weak external noise on the system dynamics, it is useful 
to introduce the notion of quasi-isolated systems [8]. A physical system is called quasi- 
isolated if it is subject to the action of very small stochastic perturbations modelling the 
random influence of surrounding. In order to explicitly show that stochastic perturbations 
are weak, they can be included in the evolution equations with a small factor a <C 1, so 
that such equations take the form 

^ y(a£,t) = v(y,a£,t) . (71) 

As a result, the multiplier matrix M(a£, t) and the Jacobian matrix J(a£,t) also contain 
this small parameter. And the local expansion exponent becomes 

X(a,t) =m|det < M(a£,t) > | . (72) 

An isolated system could be treated as the asymptotic limit of the related quasi- 
isolated system, as a — > 0. However, there is a very gentle point in taking such an 
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asymptotic limit. A system, being considered for long time, corresponds to the temporal 
limit t — > oo. The limits a — > and t — > oo may turn to be noncommuting! If these limits 
for the expansion exponent (72) do not commute, so that 

[lim, lim]X(cM)^0, (73) 

t— >oo o^u 

then the notion of an isolated system has not much sense, since infinitesimally small 
random perturbations principally change the system behaviour. Condition (73), when it 
holds true, can be accepted as a mathematical formulation for the principle of nonexistence 
of isolated systems. 

For illustration, we may employ the examples of Section 5. Thus, for the stochastic 
oscillator of subsection 5.1 we have 

X(a, t) = -Tt + a 2 ( 7l f + 7 ft 2 ) . (74) 

This gives us the limits 

lim Urn X(a,t) = — oo , lim lim X(a,t) = +oo , (75) 

where it is assumed that 72 7^ 0. For the stochastic diffusion equation of subsection 5.2, 
we get 

2 4 

X(a, t)~-?- N 3 Dt + a 2 ^— N 5 Ut + 7 2 t 2 ) , (76) 
3 5 v 7 

from where it follows that 

lim lim X(a,t) = —00 , lim lim X(a,t) = +00 . (77) 
And for the stochastic Schrodinger equation of subsection 5.3, we find 

X(a,t) = a 2 ( 7l t + 72 2 t 2 ) N, (78) 

which yields the limits 

lim lim X(a, t) — , lim lim X(a, t) = +00 , (79) 

As we see, irrespectively of whether the system without random perturbation has been 
phase conserving or contracting, it turns, in the long run, to an expanding system as soon 
as it is influenced by a noise, no matter how weak the noise is. 

In this way, at finite times the local expansion exponent can be negative or zero; for 
some systems it may, probably, fluctuate, similarly to entropy fluctuations in macroscopic 
kinetics [29]. 

But there always exist such small random perturbations that at sufficiently long times 
the expansion exponent becomes positive. This suggests to formulate the principle of 

asymptotic phase expansion, 

X(t) > (t -> 00) , (80) 
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telling that the expansion exponents of quasi-isolated systems become positive at asymp- 
totically large times. Thus, any real physical system can be treated as approximately iso- 
lated only for a finite time interval. But sooner or later, the influence of its uncontrollable 
stochastic environment will prevail and the system phase volume will start expanding. 

The principle of asymptotic expansion (80) is explicitly related to the direction of 
time. Therefore, the irreversibility of time arrow can be treated as a consequence of 
this principle. One often connects the irreversibility of time with the increase of entropy 
postulated for isolated systems by the second law of thermodynamics. The increase of 
entropy could be due to the internal chaotic nature of the microscopic dynamics [7]. 
However, as is recently reviewed by Zaslavsky [30], chaotic dynamics in real systems 
does not provide finite relaxation time to equilibrium or fast decay of fluctuations, and 
chaotic systems are not completely random in the sense originally postulated for statistical 
systems. Chaotic systems possess a property different from the regular understanding of 
randomness, a property called persistence of nonequilibrium [30]. Thus, chaotic dynamics 
of isolated systems cannot provide an explanation for the increase of entropy. Moreover, 
there is no effective definition of entropy for systems far from equilibrium [2]. Contrary 
to entropy, the local expansion exponent is defined not only for statistical systems but for 
arbitrary dynamical systems, strongly as well as weakly nonequilibrium. An important 
fact is that completely isolated systems do not exist in nature. Any real system can be 
only quasi-isolated, being subject to uncontrollable random influence of its environment. 
Although the internal chaotic dynamics, if any, may play the role, but the main cause 
for the property of asymptotic expansion (80) is the influence of stochastic surrounding. 
It is this asymptotic phase expansion that indicates the direction of time and makes 
the time arrow irreversible. The definition of the expansion exponent does not invoke 
any thermodynamic or statistical notions, but is valid for arbitrary dynamical systems. 
This is why the irreversibility of time arrow is not a privilege of only statistical and 
thermodynamic systems, but it is a common property of all evolution processes. 

7 Conclusion 

The notion of local expansion exponents is introduced being valid for arbitrary dynamical 
systems, including stochastic dynamical systems. This notion plays a fundamental role in 
the problem of pattern selection. The probabilistic approach to pattern selection yields the 
principle of minimal expansion exponent, according to which the most probable pattern 
at a given time corresponds to the minimal local expansion exponent. 

Considering the limits of large times and of weak stochastic sources, as applied to 
the expansion exponent, makes it possible to give a mathematical formulation for the 
principle of nonexistence of isolated systems, when these two limits do not commute with 
each other. 

Several examples of quasi-isolated systems suggest, as a plausible generalization, the 
principle of asymptotic phase expansion, stating that the expansion exponents of quasi- 
isolated systems become positive at sufficiently long times. 

The positive definiteness of the expansion exponent at asymptotically large times 
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implies the expansion of the system phase volume. The relation between the asymptotic 
phase expansion and time defines the direction of time and provides the foundation for 
the principle of irreversibility of time arrow, valid for all evolution processes. 
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Appendix A: Proof of Proposition 1 

When the multiplier matrix satisfies the eigenproblem (19), its elements, defined in 
equation (17), are diagonal, 

M mn (t) = 5 mn n n (t) . 
Substituting this into equation (18) gives 

Smn ^7 ^n(t) = J mn{t) l^n{t) + [/i m (t) - /jL n (t)\ <p+(t) ? 

at at 
where the property (16) has been used. For m = n, we get the equation 

^ Unit) = J nn {t)n n (t) , 

with the initial condition 

/i„(0) = 1 , 

following from condition (13). The solution of the latter equation results in the eigenvalues 
(20) containing the diagonal elements of the Jacobian matrix. In general, the latter matrix 
is not diagonal, with the nondiagonal elements being 



Jmniy) ~ 

where m ^ n. 

Appendix B: Proof of Proposition 2 

If condition (21) is valid, then equation (18) becomes 



dt 



j t M mn (t) = Y, J mk (t)M kn (t) + M mn (t) 



rf^(t) d(p m (t) 



Suppose that y? n (t) are the eigenvectors of the multiplier matrix M(t), so that the 
eigenproblem (19) holds. Then the above equation yields 

$mn(t) fJin(t) = Jmnif) fX n (t) , 

from where it follows that the Jacobian matrix is diagonal 

that is, (p n (t) are also the eigenvectors of J(t). The eigenvalues fi n (t) and J n (t) are related 
by expression (22). 

Let now (p n (t) be the eigenvectors of the Jacobian matrix J(t). Then equation (18) 
can be solved as 



M mn (t) = M mra (0)exp 



o 



dt' rmy ' dt' 
20 



dt' 



From here, with the initial condition M mn (0) = 5 mn , we have 

M mn {t) = 5 mn exp J n (t') dt'} , 

which shows that the multiplier matrix is diagonal. Hence, y?„(t) are also the eigenvectors 
of M(t), with the eigenvalues (22). 

Appendix C: Proof of Proposition 3 
Differentiating the local expansion exponent (15) gives 

j X(t) = Re Tr J(t) 

for all t > 0. Due to the existence of the limit (23), 

lim ^-X(t) = Re TrJ . 

With J n being the eigenvalues of J, one has 

Re TrJ = £ Re J n = ]T A n , 

n n 

where \ n = ReJ„ are the Lyapunov exponents. The eigenvectors of J do not depend on 
time, hence they satisfy condition (21). Then, by theorem 2, the limit of M(t) at t — > oo 
possesses the same eigenvectors as J. At asymptotically large t, because of representation 
(10) one has 

X(t)^J2^W(t)\ (t^oo). 

n 

And relation (22) gives 

In \n n (t)\ ~ Re f JJt') dt' (t -> oo) . 

Therefore 

1 

lim - In |/i n (t)| = Re J„ = A n . 

t— >oo J 

Thus, we prove equations (24) and (25). 



21 



References 

[1] D. Ruelle, J. Stat. Phys. 85 (1996) 1. 
[2] D. Ruelle, J. Stat. Phys. 95 (1999) 393 . 

[3] G.M. Zaslavsky, Chaos in Dynamical Systems, Harwood, New York, 1984. 
[4] L. Rondoni, E.G.D. Cohen, Nonlinearity 13 (2000) 1905. 

[5] J. von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton Uni- 
versity, Princeton, 1955. 

[6] V.I. Yukalov, Statistical Green's Functions, Queen's University, Kingston, 1998. 

[7] J.L. Lebowitz, Physica A 263(1999) 516. 

[8] V.I. Yukalov, Physica A 234 (1997) 725. 

[9] V.I. Yukalov, Phys. Rep. 208 (1991) 395. 

[10] V.I. Yukalov, Phys. Lett. A 284 (2001) 91. 

[11] V.I. Yukalov, Physica A 291 (2001) 255. 

[12] M.C. Cross, PC. Hohenberg, Rev. Mod. Phys. 65 (1993) 851. 

[13] V.I. Yukalov, Phys. Lett. A 278 (2000) 30. 

[14] V.I. Yukalov, S. Gluzman, Int. J. Mod. Phys. B 13 (1999) 1463. 

[15] V.I. Yukalov, Mod. Phys. Lett. B 14 (2000) 791. 

[16] V.I. Yukalov, Eur. Phys. J. B 20 (2001) 609. 

[17] C.W. Gardiner, Handbook of Stochastic Methods, Springer, Berlin, 1997. 

[18] R.L. Stratonovich, Topics in the Theory of Random Noise, Gordon and Breach, New 
York, 1963. 

[19] R.L. Stratonovich, Nonlinear Nonequilibrium Thermodynamics, Springer, Berlin, 
1992. 

[20] V.I. Yukalov, Phys. Rev. A 56 (1997) 5004. 
[21] V.I. Yukalov, Laser Phys. 7 (1997) 998. 

[22] D. ter Haar, Elements of Statistical Mechanics, Reinehart, New York, 1954. 
[23] L.D. Landau, E.M. Lifshitz, Statistical Mechanics, Pergamon, Oxford, 1958. 
[24] J.E. Mayer, M.G. Mayer, Statistical Mechanics, Wiley, New York, 1977. 

22 



[25] O. Penrose, Rep. Prog. Phys. 42 (1979) 1937. 

[26] V.I. Yukalov, Mosc. Univ. Phys. Bull. 25 (1970) 49. 

[27] V.I. Yukalov, Mosc. Univ. Phys. Bull. 26 (1971) 22. 

[28] W.H. Zurek, Prog. Theor. Phys. 89 (1993) 281. 

[29] B.V. Chirikov, O.V. Zhirov, J. Exp. Theor. Phys. 120 (2001) 214. 

[30] G.M. Zaslavsky, Physics Today, N8 (1999) 39. 



23 



